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The temporal magnetic correlations of the triangular lattice antiferromagnet NiGa 2 S 4 are exam¬ 
ined through thirteen decades (10“^^ — 1 sec) using ultra-high-resolution inelastic neutron scattering, 
muon spin relaxation, AC and nonlinear susceptibility measurements. Unlike the short-ranged spa¬ 
tial correlations, the temperature dependence of the temporal correlations show distinct anomalies. 

The spin fluctuation rate decreases precipitously upon cooling towards T* = 8.5 K, but fluctuations 
on the microsecond time scale then persist in an anomalous dynamical regime for 4 K < T < T*. 

As this time scale exceeds that of single site dynamics by six orders of magnitude, these fluctuations 
bear evidence of emergent degrees of freedom within the short-range correlated incommensurate 
state of NiGa 2 S 4 . 

PACS numbers: 75.40.Gb, 78.70.Nx, 76.75.+i, 75.30.Cr 


The two dimensional triangular lattice plays a promi¬ 
nent role in the field of frustrated magnetism. While 
it now seems clear that the conventional triangular lat¬ 
tice antiferromagnet has long-range non-collinear 120 de¬ 
gree Neel order irrespective of the spin quantum number, 
S , experiments on materials with generalized longer 
range interactions, spin-orbit interactions, as well as lat¬ 
tice distortions point to a rich variety of phases in the 
vicinity of the Neel phase 0 , 0 . In the quantum limit 
{S = 1/2) first considered by Anderson experiments 
on organic materials suggest a gapless spin liquid with a 
spinon Eermi surface may exist near the metal insulator 
transition 0 • Recent exact-diagonalization showed such 
quantum spin-liquid states can be stabilized by the pres¬ 
ence of randomness in the near-neighbor interactions Q . 
Deep in the insulating limit, neutron scattering exper¬ 
iments suggest conventional spin waves may be supple¬ 
mented or replaced by a spinon-like continuum even when 
long-range order is present @ . 

Eor larger spm quantum numbers theoretical 0 and 
experimental [ll| work points to complex slow spin dy¬ 
namics that is still poorly understood. Here we focus on 
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such dynamics in the S = 1 triangular lattice antiferro¬ 
magnet NiGa 2 S 4 H, [lH, which is distinguished by a 2D 
incommensurate critical wave vector and the potential 
for spin-nematic interactions [ 30 . Using high energy 
resolution inelastic neutron scattering, neutron spin echo, 
muon spin relaxation (/iSR), and AC magnetometry, we 
provide evidence for emergent MHz spin dynamics over 
a wide range of temperatures. We argue such slow spin 
dynamics is associated with emergent topologically pro¬ 
tected degrees of freedom inherent to incommensurate 
magnetism on the triangular lattice. 

Polycrystalline samples of NiGa 2 S 4 were synthesized 
by solid state reaction described elsewhere [l6|. Elastic- 
channel data were taken on the GPTAS and HER triple¬ 
axis spectrometers at the JRR-3 research reactor. The 
backscattering data were taken on the high flux backscat- 
tering (HEBS) instrument at the NIST Center for Neu¬ 
tron Research (NCNR) 0. Neutron spin-echo experi¬ 
ments were performed on the NG5-NSE at NCNR and 
iNSE spectrometers at the JRR-3. AC susceptibility and 
nonlinear susceptibilities were measured using a commer¬ 
cial SQUID magnetometer. Detailed descriptions of ex¬ 
periments, analysis of /iSR data, and discussion can be 
found in the Supplementary Information. 

NiGa 2 S 4 consists of neutral insulating layers held to¬ 
gether by the van der Waals force (Eig. 1(a)). This 
curtails electron hopping between layers and leads to 
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FIG. 1: (color online). Temperature dependence of (a) the 
spin correlation length ^ and (b) characteristic spin relaxation 
time T of NiGa 2 S 4 . The crystal structure of NiGa 2 S 4 is de¬ 
picted in the inset to (a). NQR data in (b) is after Ref. 0 . 
Above and below T* the muon spin relaxation rates are in 
the motionally-narrowed (oc r) and adiabatic (oc l/r) lim¬ 
its, respectively 0111 (see Supplementary Information for 
details). Vogel-Fulcher behavior from the AG susceptibility 
measurements is given by the dotted curve. The dashed curve 
in (a) and the solid curve towards T* in (b) are guides to the 
eyes. 


a strongly two-dimensional magnet where spin correla¬ 
tions do not extend beyond nearest-neighbor planes 
Indeed, a study of the effects of nonmagnetic impuri¬ 
ties indicates quantum collective phenomena in NiGa 2 S 4 
are controlled by intra-plane interactions [l8[. Within 
a layer, nickel ions with S = 1 reside on an equilateral 
triangular lattice. NiGa 2 S 4 is thus a rare realization of 
a highly two-dimensional quantum antiferromagnet. Up¬ 
holding expectations of anomalous magnetism, the ma¬ 
terial does not exhibit conventional magnetic-order down 
to temperature T = 50 mK. Neutron scattering reveals 
incommensurate spin correlations with an in-plane corre¬ 
lation length ^ that gradually increases with cooling but 
tends to 2.6 nm (7 lattice spacings) without a thermal 
anomaly [Fig. 1(a)]. 

In sharp contrast, temporal correlations feature clear 
thermal anomalies. Figure 1(b) shows the temperature 
dependence of the characteristic spin relaxation time r 
in NiGa 2 S 4 and is the central result of this work. Un¬ 



FIG. 2: (color online). Temperature dependence of the elastic 
channel data from NiGa 2 S 4 . The nominally elastic signals of 
the magnetic diffuse scattering from poly crystalline NiGa 2 S 4 
are obtained using thermal-neutron (energy-resolution 8E — 
2.96, 1.27 meV) and cold-neutron {5E — 0.37, 0.12 meV, 
50.9 /ieV) triple-axis spectrometers at the wave vector Q = 
0.58 A“^, and the backscattering spectrometer {5E — 0.8 
HeY) at Q = 0.76 [i3|. 


like normal magnets, NiGa 2 S 4 possesses an intermediate 
temperature regime with spin dynamics on time scales 
six orders of magnitude greater than the time scale set 
by the exchange constant J. Multiple experimental tech¬ 
niques underlie these data as described in the following. 
Gallium nuclear quadrupole resonanc e ( NQR) experi¬ 
ments 0 and /iSR measurements [13, Ell reveal critical 
slowing down of Ni-spin fluctuations between the Weiss 
temperature \0w\ = 80 K and a characteristic temper¬ 
ature T* = 8.5 K, where nominally elastic diffuse mag¬ 
netic neutron scattering develops. This is accompanied 
by a loss of the NQR signal between T* and 2.5(5) K in¬ 
dicating dynamic nuclear spin relaxation resulting from 


slow Ni spin fluctuations [20|, l2l|. //SR provides infor¬ 
mation for T <T* where the NQR signal is absent, and 
shows the spin relaxation rate decreases slightly below 
T* but remains finite down to 25 mK [2l|. These results 
demonstrate Ni spin dynamics persists to the lowest tem¬ 
peratures, with considerable spectral weight in the MHz 
frequency range. 


Figure 2 shows the temperature dependence of nom¬ 
inally elastic magnetic neutron scattering from poly¬ 
crystalline NiGa 2 S 4 as measured with varying spec¬ 
trometer energy resolutions. The data were obtained 
on thermal- and cold-neutron triple-axis spectrometers, 
and backscattering spectrometer The energy- 

resolution 5E sets the minimum spin relaxation time r ^ 
h/5E for spin fluctuations that contribute to intensity 
in the elastic channel. Thus the gradual appearance of 
elastic scattering below a ^F^-dependent crossover tern- 


















3 



Fourier time (ns) 

FIG. 3: (color online). Fourier time t dependence of the inter¬ 
mediate scattering function S{Q^t) / S{Q^0) of NiGa 2 S 4 . (a) 
T = 1.4 K data and (b) T = 3.4 K data were acquired on 
the NG5-NSE instrument at NGNR. (c) T = 10.2 K data are 
obtained on the G2-3-l-iNSE instrument at JRR-3. All the 
data were taken at wave vector transfer Q = 0.65 A“^. The 
dashed lines in (a) and (b) give the averages, and an expo¬ 
nential decay fit is shown as the solid line in (c). 


perature signals the development of spm correlations on a 
time scale beyond r ^ hjdE (Fig. 2) [22|. The downward 
shift in onset temperature with tighter energy-resolution 
indicates the time scale for spin fluctuations increases 
with cooling as presented in Fig. 1(b). 

The smallest value of 5E is obtained using the 
neutron spin echo (NSE) technique. The result¬ 
ing dependence of the intermediate scattering func¬ 
tion (ISF) S{Q,t)/S{Q,0) on Fourier time t, is shown 
in Fig. 3 for three representative temperatures. To esti¬ 
mate the relaxation time r, the data were fit to a sim¬ 
ple exponential decay form, exp(—t/r). For data taken 
at T = 10.2 K (Fig. 3(c)) where a slight reduction in 
S{Q^t)/S{Q^0) is visible up to t = 16 ns, the fit yields 
r ^ 0.11(5) /is with reduced = 3.971(6). A fit to con¬ 
stant ISF leads to = 4.014(6). Figures 3(a) and 3(b) 
show that below T*, S{Q,t)/S{Q,0) is constant and does 
not show any apparent relaxation. Assuming exponen¬ 
tial decay, the data place a lower limit of r = 1.3 jas 
on the corresponding relaxation time. The fact that 
S{Q,t)/S{Q,0) < 1 for early times indicates fast initial 
relaxation and the decrease of S{Q,t > 0.1 ns)/5'((5,0) 
with increasing temperature (Fig. 3) indicates a corre¬ 


sponding reduction in the fraction of spins participating 
in slow relaxation. Given the temperature dependence of 
the ISF for Fourier time ranging from 0.05 to 1.0 ns (Sup¬ 
plementary Information), we interpret the ISF as the vol¬ 
ume fraction of spins correlated in accordance with the 
critical wave vector Qc and relaxing at times beyond the 
Fourier time. The temperature dependence of the time 
averaged ISF for 0.05 < t < 1 ns is given in Fig. 4(a). 
Reflecting the behaviour of r for To < T < T*, where 
To ^ 4(1) K (Fig. 1(b)), there is a temperature regime 
where the time-averaged ISF is approximately constant 
possibly due to quantum tunneling. The time-averaged 
ISF however, increases upon cooling below To, indicating 
slow dynamics through a larger fraction of the sample, 
which eventually results in recovery of the NQR signal 
for T < 2.5(5) K [I^. 

For access to the MHz time scale we turn to /iSR. A 
description of the method employed to infer a relaxation 
time from /iSR data can be found in Supplementary In¬ 
formation. Interestingly, the relaxation times from res¬ 
onance data are considerably shorter than inferred from 
neutron scattering for T < 1.5T*. This does not however, 
represent a discrepancy since the techniques probe differ¬ 
ent aspects of the spin correlations. In conjunction, the 
results indicate critical spin-pair correlations probed by 
neutron scattering at Qc persist on a longer time scale 
than the Q-averaged correlations probed by NQR and 
/iSR. 

To investigate ultra-slow spin dynamics we use AC 
susceptibility measurements, which probe the uniform 
(5 = 0 response. Figure 4(b) gives the in-phase com¬ 
ponent x' a function of temperature for frequencies 
from 1 to 1000 Hz. While no significant frequency de¬ 
pendence is observed at the T* anomaly where the spin 
fluctuation rate is in the MHz range, frequency depen¬ 
dence is observed for T ^ 3 K where the NQR signal 
is recovered M- Below 2.5(5) K, a quasi-static inho¬ 
mogeneous internal field was inferred from the NQR ex¬ 
periments Il9l while /iSR is inhomogeneous at all tem¬ 
peratures El- These observations indicate a wide 
distribution of spin fluctuation rates, which is corrobo¬ 
rated by the distinct time scales inferred from /iSR, NQR 
and ac susceptibility data (Fig. 1(b)). That this distribu¬ 
tion is strongly (5-dependent should be expected given 
the (5-dependence of magnetic neutron scattering, and is 
witnessed by the longer relaxation times for Q = 0 AC 
susceptibility and Qc neutron scattering data compared 
to local probes (/iSR and NQR). The relatively slow fluc¬ 
tuations seen in NQR could be responsible for the quasi¬ 
static component of the //SR signal. Extrapolation of the 
Vogel-Fulcher curve (Supplementary Information) for the 
AC susceptibility shows that microsecond spin dynam¬ 
ics is present at least down to Tq. While neither the 
/iSR or NSE data exclude a truly static component, the 
/iSR data indicate significant MHz spin dynamics down 
to 25 mK [Hi. 
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FIG. 4: (color online). Temperature dependence of (a) the 
averaged intermediate scattering function *S'(Q,0.05 < t < 
1 ns)/AS'((5, 0), (t>) AC and (c) nonlinear susceptibilities, the 
coefficient of the third power (left-axis) and the fifth power 
(right-axis). The dashed curve and horizontal dotted line in 
(a) are guides to the eyes. The inset to (b) emphasizes the 
frequency dependence of x' cit low temperatures. The vertical 
dashed and dotted lines denote T* = 8.5 K and To = 4(1) K, 
respectively. 


At its freezing temperature, a canonical spin-glass is 
predicted to have diverging third and fifth-order nonlin¬ 
ear susceptibilities, xs and xs [23|. In search of possible 
spin-glass-type freezing in NiGa 2 S 4 , nonlinear suscepti¬ 
bility measurements were performed. Figures 4(c) shows 
the temperature dependence of xs and xs- Neither shows 
divergent behaviour at T* or Tq. The quadratic temper¬ 
ature dependence of the low-temperature magnetic spe¬ 
cific heat Cm 0 in NiGa 2 S 4 also contrasts with the 
conventional linear-T behaviour known for spin-glasses. 
Similar Cm ^ behaviour has been reported for 
quasi- 2 D Kagome-related frustrated magnets of deutero- 
nium jarosite 0 and SCGO [25|. The quasi-two- 
dimensional antiferromagnetic spin-correlations revealed 
by Q-dependent magnetic neutron scattering [110 dif- 
fer qualitatively from the canonical spin-glasses, where 
no Q-dependence is anticipated [23|. 

Despite the structural simplicity and considerable the¬ 
oretical efforts, no theory is so far able to fully account 
for the unusual magnetism of NiGa 2 S 4 . In the following 


we shall discuss models with potential relevance to the 
experimental findings. In Monte-Carlo simulations, clas¬ 
sical Heisenberg spins with near-neighbor interactions on 
the triangular lattice exhibit a precipitous decrease in the 
spin relaxation rate for T^/ J ^ 0.28 [l^, which may be 
compared to T*/J = 0.43 (J = 20(1) K from high-T sus¬ 
ceptibility data). This phenomenon has been associated 
with topological Z 2 -vortex ordering 0 . Slow spin dy¬ 
namics and a finite correlation length is predicted below 
Tv [i^. Aspects of this theory are in accordance with 
our findings of consecutive anomalies in the magnetic 
susceptibility and Cm, and the asymptotic behaviour of 
X 3 (Supplementary Information). However, the correla¬ 
tion length predicted in the low-T regime exceeds that 
observed in NiGa 2 S 4 by orders of magnitude. 

Lacking in this theory as a model for NiGa 2 S 4 is the 
quantum nature of Ni-spins and the incommensurate Qc- 
The semi-classical description of spin degrees of freedom 
thus cannot account for the impurity-spin size depen¬ 


dence of doped NiGa 2 S 4 lH, 13, Isl- Specifically, the low- 
temperature T^ asymptotic behaviour of Cm is preserved 
only for i nteg er impurity spins that support a quadrupole 
moment [IJ, ll5|. Indeed, a recent theory achieves slow 
spin dynamics in a near-neighbor model of NiGa 2 S 4 with 
impurity spins induced by bond-disorder with long-range 
interactions mediated by an antiferroquadrupolar bulk 
phase [HI. 

None of these models however, considers the in¬ 
commensurate nature of spin correlations in NiGa 2 S 4 
which are thought to arise from competing ferromag¬ 
netic nearest-neighbor and antiferromagnetic third near¬ 
neighbor interactions 0 . Apart from global rotations, 
the conventional 120 degree spin structure on the trian¬ 
gular lattice comes in two flavours distinguished by the 
scalar chirality of the three spins on a specific triangle. 
However, the incommensurate structure of NiGa 2 S 4 has 
three different wave-vector domains, which are related to 
each other by a Cs rotation. Monte-Carlo simulations 
show that Cs symmetry breaking occurs through a finite 
temperature first-order phase transition jsol, EH . Based 
on the inferred exchange interactions, the predicted crit¬ 
ical temperature (O. 24 J 3 ^ 6 K) is not far from T*. The 
Imry-Ma criterion [HI however, implies this transition 
does not survive under arbitrary weak random field dis¬ 
order. T* may thus be the remnant of this first-order 
transition and the anomalous MHz dynamics might arise 
from interfaces between the associated ( 73 -domains [HI • 

We have characterized megahertz dynamics in a short- 
range correlated state finding a temperature dependent 
relaxation time that is not reflected in temperature de¬ 
pendent spatial correlations. A spin state with emer¬ 
gent topological degrees of freedom that form at T* and 
freeze at lower-T may support decoupling between spa¬ 
tial and temporal two-point spin correlations. Low fre¬ 
quency (hjj J) and low temperature {ksT < J) spin 
fluctuations as reported here are not uncommon in geo- 
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metrically frustrated magnets but they are poorly under¬ 
stood. Owing to its structural simplicity and the avail¬ 
ability of high quality single crystals, NiGa 2 S 4 presents 
an excellent opportunity for progress towards under¬ 
standing slow dynamics in spin systems with short-range 
correlations. 
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